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I. INTRODUCTION AND SUMMARY 



The main cosmological constant problem (CCP No. 1 or CCPl, for short) hes in the 
apparent conflict between certain theoretical expectations and experimental facts (see, e.g., 
Ref. Ill] for an extensive review). The key theoretical expectation is that the zero-point 
energy of the quantum fields in the equilibrium vacuum state naturally produces an unsup- 
pressed effective cosmological constant A in the classical gravitational field equation. The key 
experimental fact is the observed negligible value of A. The qualifications 'unsuppressed' 
and 'negligible' refer to the known energy scales of elementary particle physics. Hence, 
CCPl motivates us to discover the mechanism which cancels the gravitational effects of this 
zero-point energy, without fine-tuning the theory. 

Several years ago, Dolgov ^ proposed a remarkable solution to CCPl by having an 
evolving massless vector field which dynamically cancels the effective cosmological constant 
A. The original Dolgov model, however, runs into two obstacles. The first obstacle [sl is 
that the steadily increasing vector field ruins the Newtonian gravitational dynamics of a 
localized matter distribution, e.g., the matter of the Solar System. The second obstacle, 
already noted by Dolgov himself, is that the expansion of the asymptotic Universe is too 
fast, inexorably diluting any standard-matter component initially present. 

Inspired by the g-theory approach [4, 5| to CCPl, two extended vector-field models have 
been constructed, which circumvent each obstacle separately . The question is whether 
or not there exists a further extended vector-field model which deals with both obstacles 
simultaneously. The present article answers this question affirmatively by doing the obvious, 
namely, by combining the two previous extended models. 

With this final extended vector-field model, we have a possible solution of the main cos- 
mological constant problem and no unwanted side effects. Indeed, the final extended vector- 
field model cancels an arbitrary (Planck- scale) cosmological constant A without fine-tuning, 
while maintaining the standard local Newtonian gravitational dynamics and providing for 
an acceptable late-universe Hubble expansion (there may still be an inflationary phase in the 
very early universe Q). But, this is only a 'possible' solution, because it is not clear if such 
massless vector fields exist in reality (having a consistent quantum theory and a mechanism 
to guarantee their masslessness) . A further caveat on this 'possible' solution is mentioned 
in Endnote [9] which is called in Sec. IV D[ 

Taking for granted that CCPl has been solved in principle, the next problem (CCP2) 
is to explain the small but nonzero value measured in the actual nonequilibrium Universe. 
This problem lies outside the scope of the present article. Some relevant remarks can be 
found in the recent review [s], which contains, moreover, a brief summary of g-theory. 
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The outline of this article is straightforward: first, the model is defined (Sec. [TTl), then, the 
homogeneous background solution is determined (Sec. and found to correspond to an 
attractor-type solution (Sec. lIVp . and, finally, the local gravitational dynamics of small-scale 
systems is shown to be Newtonian (Sec. |V]). Two appendices give mathematical proofs for 
the existence of attractor solutions in related but simpler vector-field models, namely, the 
original Dolgov model ^ and the model of our second article 3]- 

All calculations contained in this article are analytic. The main results are, first, the 
exact solution (13.81) with constants (13.91) and A-cancellation (13.101) and, second, the vanishing 
modification (I5.19ap of the weak-field gravity theory (I5.18P for a localized matter distribution 
in a perfect-equilibrium background. 

II. TWO- VECTOR-FIELD MODEL 

The model is presented in Sec. Ill Al together with appropriate cosmological Ansdtze for 
the fields. The reduced field equations are given in Sees. IIIBI and III Ci 

A. Action and Ansatze 

Consider a model of two massless vector fields, Aa{x) and Ba{x). This model is governed 
by the following effective action {h = c = 1): 

Planck 

Y R[g] + e(F4, Fb) + A + CmIq, (j)]^ , (2.1a) 

Epianck = (87^G)-l/^ (2.1b) 

with the Ricci scalar R{x) of the metric (7^/3 (x), the effective cosmological constant A, a 
generic massless matter field 0(x) with a standard Lagrange density Cm{x), and a function 
e{FA,FB) of the following two auxiliary variables Fa{x) and Fb{x): 

F^lg, A] = {Q,A? -^RA^A", Q^Alg, A] = VM„, (2.2a) 

FB[g, B] = (Q.Bf-^RB^B'^, QsbIq, Bj^V'B^, (2.2b) 

where denotes the covariant derivative (later also written as a semicolon in front of the 
relevant spacetime index). As will become clear at the very end of this article (Sec. IV Pp . 
the coupling constant G entering the reduced Planck energy fl2.1bp can be identified with 



Newton's gravitational coupling constant Gn, first measured by Cavendish. The magnitude 
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of the cosmological constant is considered to be of the order of the Planck energy, |A| ~ 
(£Jpianck)^- For the moment, we set the standard term Cm in (12 .lap to zero [possible zero- 
point-energy contributions from the field have already been included in A, assuming a 
proper (relativistic) regularization tracing back to the fundamental microscopic theory]. 
Combining the Ansdtze of our previous work , we take the following special function: 



e(F^, Fb) = (^Planck)' (^a + 6 J , (2.3) 

with numerical constants a = ±1 and b = —a. In explicit calculations later on, we will use 
the values a = —b = 1. The function (12. 3p possesses the following symmetry properties: 

which will turn out to be crucial for the preservation of standard Newtonian gravity on small 
scales (Sec. IV Pp . 



In this article, we start by considering a spatially fiat, homogeneous, and isotropic uni- 
verse. The corresponding Robertson-Walker (RW) metric in suitable spacetime coordinates 
is given by: 

(^gapixi, X2, X3, t)^ = (^diag[l, -a^{t), -a^{t), -a^(t)]) , (2.5) 

where a{t) is the scale factor as a function of cosmic time t. The usual Hubble parameter is 
defined hj H = {da/dt)/a. 

The following Ansdtze {2! for the background vector fields are consistent with the homo- 
geneous and isotropic background metric (12. 5p : 

A«(xi, X2, t) = Aoit) 5° , (2.6a) 

S„(xi, X2, X3, t) = Bo{t) 6l , (2.6b) 
involving only two functions of t. 

B. Reduced vector-field equations 

The variational principle for the vector fields of action (12. lap gives the following two 
equations: 

(V„Q3A + 2 ^ "^"j 9^ + V„ ( ^ j = , (2.7a) 
(V,Q3B+ 2^^-) ^ + Q3bVJ =0. (2.7b) 
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The Ansdtze (12. 5p and (12. 6 p reduce the eight partial differential equations (I2.7ap and 
f l2.7&P to the following two ordinary differential equations (ODEs): 

(i„ + 3ffi.-6ffM,)^+(^, + 3i/^„)^(^) =0, (2.8a) 

B„H-3i/B„-6i/-^B„)^H-(B„H-3i/B„)^(^) =0, (2.8b) 

where the overdot stands for differentiation with respect to the cosmic time t. Similarly, the 
g-theory-type variables from (12 ■2p become 

QsaW = j^Aoit) + 3HAoit), (2.9a) 

Q3B{t) = j^Boit) + 3HBo{t). (2.9b) 

Observe that having Aq oc t and H oc 1/t makes Q^a in (I2.9ap into a genuine (spacetime- 
independent) g-theory variable j4| and similarly for Q-^b in (12.9bp . 



C. Generalized FRW equations 



The energy-momentum tensor of the vector fields is calculated by varying the action with 
respect to the metric tensor Qais- This tensor is found to be given by 



T, 



a/3 



e{FA,FB)-2FA 



de 
dTA 



-2F, 



B 



de 



OFa 

-v„v^ 



Rc^pA^ -RA^Ap] + 



A^ TT^ + 



de 
'dF^ 
de 



9a(3 



Rap B^ — R BaBjS 



dFB 



2 1 ,2 de de 



(2.10) 



dFA dFBj V dFA dFB, 

At this time, we also introduce a contribution to the total energy-momentum tensor from 
the standard-matter sector of the theory, that is, we consider having Cm 7^ in the origi- 
nal effective action (12. lap . In the cosmological context, the standard-matter component is 
described by a homogenous relativistic fluid. Note that this physical setup is not altogether 
unrealistic, as the masses of the standard-model particles are negligible for temperatures 
T ~ Epianck > GeV. 

From the previous Ansdtze (12. 3p . (12. 5p . and (12. 6p . the generalized Friedmann-Robertson- 
Walker (FRW) equations and the standard-matter energy-conservation equation are 



3H' 



2H + 3H^ 



PM 



{F 
{F 



Planck J 



\-2 



Planck J 



A + p{A, B) + pm 



A-P{A,B)-wmPm 



-3{1 + wm)H Pm , 



(2.11a) 

(2.11b) 
(2.11c) 
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where the last equation describes the adiabatic evolution of a perfect relativistic fluid with 
a homogeneous energy density pM(t) and pressure PAiit) = wm Puif) for constant equation- 
of-state parameter wm = 1/3. The vector-field energy density (from Tq^ = p) and isotropic 
pressure (from Tj * = —PSj'^) appearing in (12. lip are given by 

p{A, B) = e{FA, Fs) + 3{H + 3 H') ^ + B', ^) 

P{A, B) = -e{FA, Fs) + {H + 3 H^) {aI ^ + Bl ^) 

d f d e „9(9e\ d^ f d e d e \ /^.^-.x 

''Jt + ^» Mi ) - a?J + ^» aTi j • 

where the symmetry property (I2.4ap has been taken into account. As (I2.12ap contains a 
term if, for example, it is clear that (12.1 lap is not the standard Friedmann equation. 



III. ASYMPTOTIC SOLUTION 



It is a straightforward exercise to determine the asymptotic {t — )■ oo) solution from the 
reduced field equations as given in Sec. [Tll In a first reading, it is possible to skip the 
technical details and to jump ahead to Sec. IIII C\ which contains the main physics result of 
this section. 



A. Dimensionless ODEs 

As in our previous articles 0, 0] , we introduce dimensionless variables by rescaling with 
appropriate powers of the reduced Planck energy -Epianck without additional numerical fac- 
tors. Specifically, we replace 

{A, e, t, H] {A, e, r, h] , (3.1a) 

{QsA, QzB, ^0, ^0, Pm] {gsA, gsB, V, w, Tm} ■ (3.1b) 

The following dimensionless ODEs for the vector fields v{t) and w{t) and the Hubble 
parameter /i(r) result from the previous vector-field, generalized Friedmann, and matter 



6 



energy-conservation equations: 

= (fi + 3fti-6ft^,)|-+(i + 3;..)|:(|-) , (3.2a) 

0= (* + 3ft^-6;.^u-)^+(^ + 3;..)A(^). (3.2b) 

= 3 /i^ - A - e - - 3 ^/i + 3 /i^j qab - 3 /i ^ qab , (3.2c) 

= rM + 4:hrM, (3.2d) 

where the overdot now stands for differentiation with respect to the dimensionless cosmic 
time r. In addition, we have the following definitions: 

e = Ia/ Ib - Ib/ Ia, (3.3a) 

e = e-2U [del Of a) - 2 Jb [de/OfB] = e , (3.3b) 

= (v + 3hvY + 3{h + 2h^)v'^ , (3.3c) 

= {w + 3hwy + 3{h + 2h'^)w^, (3.3d) 

gAB = [de/dfA) + (de/dfB) . (3.3e) 



A further FRW equation, given by the dimensionless version of fl2.11bp . can be shown to be 
consistent with the above ODEs. 

Using the symmetry property (I2.4al) in fl3.2al) and fl3.2bl) . it can be shown that bound- 
ary conditions at r = tq with v{to)/w{to) = v{tq) /w{to) give proportional v{t) and w{t) 
solutions: v{t) = [v{tq) /w{To)]w{r). 

B. Expansion coefficients 

The asymptotic solution of the differential Eqs. (13. 2p . for A of arbitrary sign, is given by 
the following series: 

v(t) = ao T + ai + a2 + 0(r~^) , (3.4a) 

w{r) = f3oT + f3i + f32 T-^ + 0(r-2) , (3.4b) 

h{r) = 70 + 71 r'^ + 72 r'' + 0(r-^) , (3.4c) 

rM(r) = 60 + 6i r'^ + 62 r'^ + OiT~') , (3.4d) 
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with leading-order coefficients: 

cto = 1 , (3.5a) 

/3o = ±7(A/2) + v/l + (A/2)2, (3.5b) 

70 = 1/2, (3.5c) 
5o = 3/4, (3.5d) 

next-to- leading-order coefficients: 

«! = — 2ao7i, (3.6a) 

(3i = -2/3o7i, (3.6b) 

71 = 7i, (3.6c) 
^1 = 3 7i , (3.6d) 

and next-to-next-to-leading-order coefficients: 

02 = , (3.7a) 

P2 = 0, (3.7b) 

72 = 2 (71)' , (3.7c) 
52 = 9{liY. (3.7d) 



These vector and metric fields have only one arbitrary constant, 71, which we interpret 
as being due to the time-shift invariance of the equations (r — > r + const). The general 
(attractor-type) solution of the three second-order ODEs and the single first-order ODE in 
(13. 2 p will have seven arbitrary constants (see Sec. IIVD|) . 

Different starting values of f (t), w{t), hij), and ruij), at large enough r = Xstart and in 
an appropriate domain, give different values of 71. Excluded starting values are those with 

{f (Tstart), t'(rstart)} = {0, 0} aud/or {w(rstart), w(rstart)} = {0, 0} aud/or rAf(rstart) = 0. 

C. Dynamic cancellation of A 

The calculational details of this section and the next should not make us forget that the 
vector fields of the model cancel the effective cosmological constant A exactly and without 
fine-tuning. 
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Indeed, the field equations (12. 8 p and (12. lip give nonzero vector-field components, a Hubble 
parameter, and a matter energy density of the form (13. 4p for coefficient 71 = 0, 

Ao{t) = ao (^pianck)^ ^ , (3.8a) 

Bo{t) = /3o (^Planck)' (3.8b) 

Hit) = 70 , (3.8c) 

PM{t) = 60 (^Planck)' t"' , (3.8d) 

where the overall normalization of Aq and Bq is irrelevant, as only the ratio of the vector-field 
components enters the action (12. lap according to Eqs. (12. 2p . (12. 3p . and (12. 6p . The coefficients 
Q^O; /3o) 70) and 60 in (13. 8 p are not put in by hand but appear dynamically. Specifically, the 
following values have been calculated in Sec. IIIIBI 



(/3o/«o)' = I A/(^Planck)' + y 1 + J AV(^Planck)« , (3.9a) 

70 = 1/2, (3.9b) 
^0 = 3/4. (3.9c) 

These particular fields give an exact cancellation of A appearing on the right-hand side 
of the generalized FRW Eqs. fimi) . 



A + e(FA, Fb) = a + e({QsA)^ {Q3 

\ / equil V 



. , - - , , -3BJ 
equil V 

= A + (Epianck)' [(«o//3o)' - (/3o/«o)'] = , (3.10) 

where the definitions (12. 2p . (12. 3p . and (12. 9p have been used for Fa,b, e, and Q3a,b, respec- 
tively. With the nullification f l310|) . the FRW Eqs. (12TT]) are solved to order t^. 

Including the higher-order terms of the asymptotic solution (13. 4p . also called the perfect- 
equilibrium solution later on, we can evaluate the effective vacuum energy density of what 
may be called the microscopic dark-energy component, that is, the energy density not from 
standard matter but from the initial ('bare') cosmological constant, the vector fields, and 
the modified gravity. A convenient definition for a spatially fiat RW universe is as follows: 

PV'-micro(^) = 3 (-Epianck)^ H {1^ — Puit) , (3-11) 

which was simply denoted pv in Ref . 0] • The result from the asymptotic solution (13. 4p is 
py-micro(t) = 0(t-^), (3.12a) 

asymp. sol. 
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which imphes 

hm yOy-micro {t)/pM{t)\ =0. (3.12b) 

i— >c» asymp. sol. 

Result ( I3.12bl) traces back to the special properties of the e-function (12 .Sp and was absent for 
the simpler models of Refs. jol, 7|, which exhibited the behavior pv'-micro(^) oc t^^. Assuming 
the relevance of our model function (12. 3 p to physics, the implication is that a new mechanism 
is needed to explain the observed finite remnant vacuum energy density of order (meV)'^. 

Expanding on the last remarks of the previous paragraph, it is not difficult to see what 
the implications are for the present energy-density ratio of dark energy and matter. For the 
sake of the argument, use |pv'-micro(^)| = which may still be an overestimate as quantum- 
dissipative effects can be expected to produce an exponential decrease (cf. Ref. jsj and paper 
[15] quoted therein). A present cold-dark-matter energy density of the order of the critical 
energy density gives pcDM(to) ~ ^piLck for ^Planck = 1/-Epianck ~ 10"''^ s and to ~ 
{cHq)^^ ~ 10^'^ s. The present energy-density ratio would then be completely negligible, 
|py-micro(^o)|/PcDM(^o) ~ (^pianckAo)^ ~ 10"^"^*^. In fact, the ratio would already be extremely 
small near the electroweak crossover: |pv'-micro(tew)|/pM(tew) ~ (-£^ew/-Epianck)'' ~ 10~^°, for 
E 

cw ~ TeV and tgw ~ -^Pianck/(-^ew)^) derived from the spatially-flat Friedmann equation 
with pm ~ (-Eew)^- With negligible pv-micro(t) from the microscopic variables (A^, i?^, and 
effectively A), further contributions to the vacuum energy density pv'-macro(t) may come from 
phase transitions and mass effects of the macroscopic standard-model fields. As discussed 
in Ref. jst, the resulting pv-macro(i':) may decrease stepwise, approximately as tpianck 

In conclusion, the exact solution (13. 8p is of paramount importance, especially if it is an 
attractor-type solution. This attractor-type behavior will be discussed in the next section. 

IV. ATTRACTOR-TYPE SOLUTIONS 

The present section is a direct follow-up of the previous one and is also rather technical. 
In order to get an idea of the attractor-type behavior, it is possible, in a first reading, to 
consider only Sec. IIVCI 

A. Mathematical considerations 

The model of interest has an action-density term e(F^,FB) as given by (12. 3p . For com- 
pleteness, two simpler models are discussed in the appendices: in App. A, the original Dolgoy 
model with just a {QiY term in the action density and, in App. B, our previous model 
with a single Fa term as defined in (I2.2ap . 
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It turns out, however, that the first-order system of differential equations for the ^{Fai Fb) 
model does not have the relatively simple structure as found in the appendices, specifically, 
Eqs. (]A12aP and (]Blla|) . Physically, the extra complications may be due to the fact the 
e{FAi Fb) model is really an f{R) modified-gravity theory, which entails higher-derivative 
field equations [in our case, f l2.11b|) has third-order derivatives of A^it), Boit), and H(t)]. 

One possible way forward would be to rewrite this particular modified-gravity theory 
as a scalar-tensor theory (more precisely, a scalar- vector-tensor theory). Instead, we prefer 
to adopt a low-tech (read brute-force) approach by pushing the explicit solutions as far as 
possible. This approach suffices to show the attractor-type behavior, even though it lacks 
mathematical rigor compared to the approach in the appendices. In fact, what would be 
needed here is the mathematical proof that the infinite sums in the expressions of Sec. IIVDI 
converge, but we will simply assume this to be the case, as has been done in most of the 
literature on the subject (cf. Ref. 0, b]). Still, awaiting this rigorous proof and the precise 
knowledge of the attractor domain, we will only speak about 'attractor-type solutions' of 
the e{FA, Fb) model rather than 'the attractor solution' tout court. 

B. ODEs 

The complete system of differential equations from Sec. IIII Al can be written as follows: 

de 



= V + 3 hi) — Q h"^ V + (v + 3 hv) — ha. 

^ ' dr 



= iv + 3 hw — Qh? w + iw + 3hw) -^hi 

^ dr 



dfA 
de 



(4.1a) 
(4.1b) 



dfB 

= 3 h'^ - X - e - tm ~ 3 {h + 3 h^) gAB ~ 3 h QAB , (4.1c) 

= 2h + 3h^ - X-e + ^rM+ {h + 3h^)gAB-2hgAB-gAB, (4.1d) 

= rM + 4:hrM, (4.1e) 

where e, /a, /b, and gAB have already been defined in (13. 3p . We will now give several 
explicit analytic solutions of these ODEs. 
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C. Particular class of exact solutions 



The differential system (14. ip has the following class of exact solutions for r > tq : 

vir) = (r - To) Ci + , (4.2a) 

w{t) = (r - To) C, + Q , (4.2b) 

Ci (r - ro)^/^ 

h{r) = -i/^ , (4.2c) 
3/4 

rM(r) = ' , (4.2d) 

(^ - TqY 

with a real constant tq G M, nonvanishing real constants Ci, C2 G M\{0}, a real constant 
C3 G M, and a discrete constant C4 G { — 1, +!}• The constants tq, C3, and C4 are arbitrary. 
The real ratio C1/C2 is determined by the input cosmological constant A via a quartic 
equation, 

C1/C2 = Rc, (4.3a) 

(i?c)' + A(/?c)' = 1, (4.3b) 

as follows from, e.g., the generalized Friedmann equation f l4.1cp . Hence, the number of free 
parameters in (14. 2 p is four: tq, (Ci C2), C3, and C4. The physically relevant parameters are, 
however, only the ratio Cs/i^Ci C2) and the relative sign C4. 

Observe that all solutions in f l4.2p give for the effective vacuum energy density of the 
microscopic degrees of freedom an exactly vanishing result, 

Ci/C2=Rc 

rv-micro(r) = 0, (4.4) 

TO: C*!, C2, Ci, C4 

with definition rv/-micro(T) = 3 hijY — TMij) from (13. lip and Rc the positive or negative real 
solution of ( ]4.3bp . Result (14. 4 p also holds for the special case C3 = 0, which corresponds to 
the perfect-equilibrium solution (13. 8p with constants (13. 9p and an arbitrary time-shift. 

For the case of C3 7^ and C4 = —1, the rescaled solutions v{t)/Ci and w{t)/C2 in (14. 2 p 
are different at finite values of r, specifically, v{t)/Ci — w{t)/C2 oc (r — tq)^^/^. Still, both 
of these functions v{t)/Ci and w{t)/C2 approach the same asymptotic solution, the one 
from above, the other from below. This is precisely the attractor-type behavior discussed 
in Sec. IIV Al and the two appendices (see also Refs. 0, (c)] and jsl, Q] for related numerical 
results) . 
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D. Series and attractor-type behavior 



A generalized Ansatz for a nontrivial solution of (14. ip at r > ri > is as follows: 

-1 , , h/n] + (r - (r - n) [v,/n - v^/rf] + (r - ri)2 

T v{t) = -'^ 



OO / 

n=l ^ 



r3 



2\ " 



(4.5a) 



(/3or)-^«;(r) 



hi/(/^o^i)] + (^-^i)' 
l + (r-ri)2 

, (r - n) [w^/{l3o ri) - «;i/(/3o rf)] + (r - n)^ 



n.=l ^ 



1 + (r - ri)3 



r3 



(4.5b) 



2r /i(r) 



[2ri/ii] + (r-ri)2 (r - n) [2 n /^i + 2 /ii] + (r - n)^ 



r - ri)2 1 + (r - ri)3 

r-Ti)^^ 



(4/3) rM(r) 



n=l ^ ^ 

[(4/3)rfrMi]+(r-n)^ 
1 + (r - ri)2 
(r - n) [(8/3) (1 - 2 n/ii) n Tmi] + (r - 



OO ^ 

E (- 

n=l ^ 



T - Ti 

r3 



1 + (r-ri)3 

\2~ 



(4.5d) 



with /3o given by ( 13. 5p . The seven constant parameters fi, wi, wi, /ii, /ii, and r^i in 
(14. 5p represent the initial values of the functions and their first derivatives at r = ri: 

v{ti) = vi , v{ti) = vi , (4.6a) 

w{ti) = wi , w{ti) = wi , (4.6b) 

h{n) = h, , h{n) = h, , (4.6c) 

rMin) = tmi, (4.6d) 

where rM{T) requires only a single boundary condition value as its ODE is first-order, the 
other ODEs being second-order. These initial values must be sufficiently close to those of 
the perfect-equilibrium solution, given by ( 14. 2 p with Ci = 1 and C3 = tq = 0. 
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Inserting the expansions (14.51) into (14. ip gives values for the coefficients a„, 6„, c^, and 
dn in terms of the initial conditions Vi, . . . , vmi- The expressions for these coefficients are 
rather bulky (even for r ^ ri) and, here, we only indicate the dependence on the initial 
conditions for the ffist few coefficients, 

(4.7a) 
(4.7b) 
(4.7c) 
(4.7d) 
(4.7e) 
(4.7f) 
(4.7g) 

In the limit of large cosmic times (that is, large on the scale of the Planck time, r ^ ri), 
the corresponding solution takes the following form: 

v{t) = r - 1 - ci + 0{t-^) , (4.8a) 

w{t) = /3o (r - 1 - ci) + 0(r-5) , (4.8b) 



ai 


= ai(ri, hi, 


hi), 








bi 


= 6i(ri, hi. 


hi). 








Ci 


= Ci(ri, hi. 


hi). 


for i - 


= 1, .. 


.,5, 


di 


= di{Ti, hi. 


hi). 








aj 


= cijiri, Vi, 


Vi, hi, hi) , 


for j 


= 2, . 


..,7, 


6, 


= bj{Ti, Wi, 


Wi, hi, hi) , 


for j 


= 2, . 


..,7, 


dk 


= dkin, hi 


hi,rMi) , 


for k 


= 2, . 


..,5 



h{r) = — 



i±£l)! + of4 



rnir) 



3 
4^2 



n=0 
4 

^(n + 1) 

n=0 



:i+ci)" 



T" 



Extrapolating this result, we obtain the asymptotic (perfect-equilibrium) solution, 

T-Ti, 



'^asymp(''"j 



W 



asymp V 



^asymp(^) 



asymp 



(4.8c) 
(4.8d) 

(4.9a) 
(4.9b) 
(4.9c) 

(4.9d) 



where ti = 1 + ci. Observe that, apart from the overall time-shift ti, the obtained asymptotic 
solution is independent of the initial conditions (14. 6 p encoded in the Ansatz (14. 5p . 

The tentative conclusion is that different initial conditions give different solutions, which, 
however, approach the same asymptotic solution ( 14. 9p . Hence, there is an attractor-type 
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behavior. But, as explained in Sec. IIV A| this conclusion needs to be proven rigorously and 
the proper attractor domain needs to be determined. 

V. SECOND-ORDER PERTURBATIONS 

We, now, turn to localized perturbations of the metric tensor field and the two vector 
fields. Denoting the four spacetime coordinates (xi, X2, X3, t) collectively as x, we consider 
the tensor field 

gadix) = gal3{t) + hal3{x), (5.1) 

with the metric Qapit) of the fiat RW spacetime (12.51) and hap{x) the perturbation 
(l^a/?! 1)- On small scales, the relevant background metric is the standard Minkowski 
metric rjaii = diag(l, —1, —1, —1). In addition, we consider the two vector fields 

Aa{x) = Ao,{t) + 5Aa{x) , (5.2a) 

Baix) = Ba{t) + SBaix) , (5.2b) 

with A^{t) = Ao{t) 5°, B^{t) = Bo{t) (5°, \SA^\ < \Ao\, and \6B^\ < \Bo\. 

This section is highly technical and, in a first reading, it is possible to skip ahead to 
Sec. |VD] with the main physics result of this section. 



A. Variation of the vector-field Lagrange density 

The Lagrange density of the vector fields is given by Ca,b = A + ^{Fa^Fb), where the 
effective cosmological constant A has been included for convenience. To second order, the 
perturbed Lagrange density reads 

/-(perturb.) _ ^(0) ^(1) ^(2) /p. 

'~'A,B — '~'A,B "T '~'A,B "T '~'A,B ' W-"^) 

with 

4% = A + e{FA,FB), (5.4a) 

r!;i^^.<"F. + ^^<"f.. (5.4b) 
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The first- and second-order variations of Qa and Fa are 

S^'^QsA = 5 A': + h'^P A^,p - g'^^d^^'^VlpA^ , (5.5a) 
S^'^QsA = h^^ SA^;p - g'^^S^'^Tl^SA, - h^^ S^'^Tl^A, - g^^ 5^^^TlpA, , (5.5b) 

= 2QsaS'^'^Q,a - ^i?(2A"M, + h'^^ A^Afi^ - ^A'S^'^R , (5.5c) 
5^^^ Fa = 2Q,A5^'^QsA + iS^'^Q3Ay-lR[5A''5A^ + 2h''^AJA^'^ 

-^5(^)i?(^2A"M„ + h''^ Aa^Ap^ - ^A^S^^'^R . (5.5d) 

Replacing A^ and 6Aa in (15. 5 p by and 6Ba gives tlie first- and second-order variations 
of Qb and Fb- 

For future use, we rewrite C^}^^ and C^ab dimensionless form, 

4S = QIao ^f'fA + QIbo j^ffB , (5.6a) 

where S^^^ fA,B and 6^'^^fA,B correspond to 6^^^Fa,b and 6^'^^Fa,b expressed in terms of dimen- 
sionless variables Va, Wa and Va, Wa- These dimensionless variables are defined as follows: 

v^{t) = A^{t) , W,(t) = B^{t) , (5.7a) 

V3A0 

^^a(a;) = 7r^MQ(x), m;c,(x) = — ^(55c,(x) , (5.7b) 

2a = Va-Wa, (5.7c) 

with dimensional constants Qsao and Qs^o- In fl5.7cp . we have added the definition of 
2q,, which will be used extensively in the next subsections. Note also that, in the above 
definitions, the background fields are distinguished by a bar and the perturbation fields by 
a hat. 
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B. Equations for the vector-field perturbations 



The equations of the vector-field perturbations are 



+ 



1 dt 



2 OFa 



2 2 
SAO ^fA + QiBO 



OFaOFb 



OFa 

SfB = 0, (5.8a) 



2 2 t^^e 



OFaOFb 



9^ r 



dFr, 



+ 



1 <9e 



2 



B 



3B0 Qp2 + QIao 



oFAOFB^^^r^-'^ ^^-^^^ 

with 5fA,B = S^^^fA,B, Sq3A,3B = ^^^^qsASB^ 5R = 5^^'^ R. Furthermore, we have q^A = 
tJo + S/iIJo and q^^ = Wq + 3hwo. Note that the above equations for the perturbations Va 
and Wa carry a third derivative of the metric perturbation hajs, since 5/^ and 5/^ contain 
(5-R, which aheady has a second derivative of hai3- 

Using the background vector-field Eqs. fl2.7ap and ( \2.7b\i . the perturbation Eqs. f l5.8ap 
and ( 15.8&|) can be reduced to 



da 1 SVIa + 
da{6QB + 



5q 



'3A 



q^A 

5q3B 

qsB 



+ 



1 



+ 



2g3A 
1 



/3B 



RVa + 
RWa + 



6R-R 
6R-R 



5q 



'3A 



93A 

^g3B 
^3B 



0, 



with definitions 



B 



de \ 


' fd^e 


SFa + 


d'^e 


SFa) 


[dFl 


OFaOFb 


de \ 


-Wd^e 


6Fb + 


d^e 


OFb) 


[dFl 


OFaOFb 



SF 



B 



SF, 



(5.9a) 
(5.9b) 

(5.10a) 
(5.10b) 



Taking f ^(t) = Wa{t) = (ait) for (ait) = (C(t), 0, 0, 0), we find q^^ = q-^^ = q?, = C+3 h (. 
Subtracting fl5.9&P from ( I5.9ap then gives 



2q,V%-RC^zp +-Vf'zp) + —Rza-——.RCaV% = 0, (5.11) 



^3 

where z'a has already been defined in (15. 7p and 



2^3 



2 (^3)^ 



{qs? 



{Q. 



3B) 



de 
Wa 

de 



OF, 



B 



dFl 
dFl 



de 



-1 



d^e 



dFsJ dFAdFs 



de 



d^e 



dFA dFAdFB 



(5.12) 
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Notice that in f l5.1ip is not coupled to the metric perturbation h^p : z^ depends only on 
the functions C,{t) and H{t) from the background fields, together with the initial conditions 
for Va and Wa- This result follows from the symmetry properties of the function e(Fyi, -Fg). 
Substituting the e function (I2.3p into fl5.16p . we find 

„ _ 1 

In the perfect-equilibrium state with Hubble parameter H{t) = 1/2 t^^ (implying R = 
0) and constant values of Qsa and Qsb (as mentioned in the last sentence of Sec. Ill Bp . 
Eqs. (15.1ip and ( I5.13P give the following final equation for = Va — Wa- 

Wa^^zp = dalt-'^^d^e'^zp)] =0. (5.14) 

equil. background 

For perturbation fields which are analytic and of finite support {v^ = = for t G [0, T] 
and \x\ > R), the solution is trivial and 

local perturb. 

Za =0. (5.15) 

cquil. background 

In other words, the two linear vector- field perturbations turn out to be equal, 6Aa{x) = 
6Ba{x), which is the same result as obtained in Ref. gI by different methods. The explana- 
tion of fl5.15p is simple: the localized perturbation fields Va and Wa obey the same equation 
and their boundary conditions over an exterior region are also the same (zero, in fact). 

C. Energy-momentum tensor of the vector-field perturbations 

The linear perturbation of the energy-momentum tensor of the vector fields is given by the 
following expression (only the arguments h, v, and w are shown explicitly on the left-hand 
side): 

Qafi[h, V, w] = {A + e) hap 

+ 1 (/iAO - /iBo) ( (2 93 V^Sa - ^ Zx) gap + 2 Rap C^Zx-R {Ca Zp + Cp Za) ) 

+ ( V„ - gap (2/iBo zx + ubo (2 ^3 W% - RC^ zx) C') 

mi^Ao-i^Bo) (2q^\/^zx-Rezx^ (^RapC-RCaCp), (5-16) 
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with definitions 



f^A = (QsaY TTTT, f^B = {QsbY TTTT , (5.17a) 
^ (Q:,Af [(Q,a)' + (Q:,Bf 9^^) . (5.17b) 

^.^(Q.s?{(Q.sY^^HQ..fJ^), (5.17c) 

SO that hao + /Ubo = and z/^o + t-'bq = for the special function ( 12. 3p and the perfect- 
equihbrium background fields with a Ricci-fiat spacetime {R = 0). 

Manifestly, Qai3[h, v, w] does not contain derivatives of the metric perturbation hajs 
and depends only on the difference between the vector-field perturbations, 2q, as defined 
by f l5.7cp . These results rely on the symmetry properties f l2.4ap and f l2.4&l) of the special 
€{Fa, Fb) function (12. 3p and on the fact that Fa and Fb are quadratic with respect to the 
vector fields and that the background fields evolve identically as mentioned below fIS.lObp . 



D. Standard local Newtonian dynamics 

With the results of the previous two subsections, we can, at last, turn to the physical 
question of interest: the gravitational self- interaction of small (noncosmo logical) systems. 
This has been discussed extensively in our previous article [g], so we can be brief. 

The linear equation for the weak gravitational field from a localized matter distribution 
is then 

UKp + lQTiGSo,p = Q, (5.18a) 
So.p = T^p-\vapV^'T^5, (5.18b) 
= Tjp^n , (5.18c) 

where the harmonic gauge, dah'^^ = (1/2) 9/? , has been used to simplify the standard 
derivative term on the left-hand side of (I5.18ap . with d'Alembertian □ = 7]""^ dad p. The 
only new contribution appears as the second term on the right-hand side of (I5.18cp and has 
been given in (I5.16p . 

Several comments are in order. First, note that the background fields Vai'T) and ^^(r) 
are such that the A + e term in f l5.16p vanishes for the perfect-equilibrium background; see, 
in particular, the derivation f l3.10p . Second, recall that the energy- momentum tensor 
of the perturbations depends only on the metric perturbation hap (but not its derivatives) 
and the difference of the vector-field perturbations. Specifically, the behavior is as follows. 
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in a symbolic notation: Qai3[h, v, = Qai3[h, (v — w), (V + + V^) {v — w)^ The 
main input for this result is that the normalized background vector fields evolve identically, 
Va = Wa = Ca(t) for t — )• oo. But this is precisely what was found in Sec. Illli The evolution 
of = Va — Wa is, therefore, not affected by the metric perturbation hai3 (at least, to the 
linear order in perturbation theory considered). Moreover, f lS.lSp states that Sq, vanishes 
due to the boundary conditions at infinity (the energy density of the matter perturbation 
being localized in space and time). 

With A + e = and 2^ = nullifying (15.161) . the conclusion is that the nonstandard term 
in (I5.18P drops out, 

^ -| local perturb. 

Qai3[h,v,w] =0, (5.19a) 

J equil. background 

and that the linear weak-gravity field equation (in harmonic gauge) equals the one of general 
relativity, 

^ / ^ 1 local perturb. 

Dha^ + WirG 5^^^""'^ = , (5. 19b) 

J equil. background 

with the standard-matter source term 

Q (matter) _ ryn (matter) 1 yS rp {ma.tter) {KTQn\ 

•^al3 = -'a/J - -^Val3V -y6 ■ I^O.iyCJ 

As mentioned before, these results hold for perfect-equilibrium background fields [given by 
(12. 5p , (12. 6p , and (13. 8 p in dimensional form or (14. 9 p in dimensionless form] , which have dynam- 
ically canceled the cosmological constant A (see Sec lIIICj) . Recall that the main cosmological 
constant problem, CCPl as formulated in Sec. [H is precisely concerned with the dynamic 
cancellation of A in the equilibrium state of the quantum vacuum. The study of small 
self-gravitating systems in a nonequilibrium background (even if this background rapidly 
approaches the equilibrium state, as discussed in the penultimate paragraph of Sec. IIII CO . 
lies outside the scope of the present article j^. 

Equation (]5.19bp shows, in particular, that the standard Newtonian law of gravity (i.e., 
the Poisson equation) holds for local nonrelativistic matter distributions such as the Solar 
System or the Galaxy. This implies that the constant G in (]5.19bp . which traces back to 
the original action (12. ip . can be identified with Newton's gravitational coupling constant, 

G = Gn = 6.6743(7) x 10~^^ m^ kg"^ s"^ , (5.20) 

where the numerical value has been taken from the CODATA-2006 compilation 13|. The 
corresponding numerical value of the gravitational energy scale defined in fl2.1bp is then the 
usual one, -Epianck ~ 2.44 x 10^^ GeV. 

The present article, just as its predecessor j6|, only considers the linear theory of small 
self-gravitating systems. It remains to be seen whether or not the present setup reproduces 
locally the standard nonlinear theory, i.e., general relativity. 
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Appendix A: Attractor solution in a model with a {Qi)^ term 
1. ODEs and new variables 



The original Dolgov model [2[, with a single massless vector field Aa{x) and a positive 
cosmological constant A, is defined by the action (12 .lap , setting Ba{x) = and having a 
vacuum-energy-density term eo based on a different contraction of the vector-field deriva- 
tives, 

eD = iQi)\ (Ala) 
iQ,r ^ A^,p A-^^ = \3H'Al (Alb) 



where the last step in (lAlbp holds for the RW metric (12. 5p and the isotropic Ansatz (I2.6ap . 



With the dimensionless variables (13. ip . the basic equations are given by the reduced 
vector-field and FRW equations: 

V + 3hv -Sh'^v = 0, (A2a) 

3/^2 = X-ij^-Sh^y\ (A2b) 

2h + 3h^ = X-v^ -Sh'^v^ -2hv^ -4:hvv + 2v^ , (A2c) 

for A > 0. It is possible to set A = 1 by an appropriate rescahng of r and 1/h, but we prefer 
to keep A explicit, in order to facilitate comparison with the models of App. [B]and Sec. IIV[ 
The system of differential equations can now be rewritten as follows: 

V + 3hv -3h'^v = 0, (A3a) 

3{l + v^)h^ - \ + v^ = , (A3b) 

{l + v'^)h-v'^ + 2hvv = 0. (A3c) 

Next, introduce new variables (using the natural logarithm 'In'): 

yi = V , 1/2 = hv , (A4a) 

s = ln(a) - In(astart) , (A4b) 
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where a = a(r) is the scale factor of the flat RW metric and Ostart its value at r = Xstart- 
Writing flA3|) in terms of the new variables (IA4|) gives the following autonomous system of 
differential equations: 

y[ = Fiivi, y^) , y'^ = F^ivi, y2) , (A5) 
where the prime stands for differentiation with respect to s and 

Fi{yuy2) = -3(?/i-?/2), (A6a) 
F2{yi, 2/2) = (a - + 3 2/2 - 6 yl) . (A6b) 

Recall that the system (lASP is called autonomous because the independent variable s does 
not occur explicitly (see, e.g., Refs. 0, Q for background material). 

2. Critical points 

A critical point (?/io, 2/20) of system flASp is defined as follows: 

Fi (1/1,1/2) =i^2(yi, 1/2) =0. (A7) 



j/10,3/20 



yio, j/20 



A straightforward calculation gives two such critical points, 

2/±=±v^/2, z/2l) = ±v^/2, (A8) 
corresponding to the asymptotic solutions 

^alymp = ± (v^ / 2) T, /l^^ymp = ^"^ , (A9) 

in terms of the original variables.^ 



3. Stability analysis: Linearization 

Make the following shift of variables: 

Vx = yio + Yi, y2 = 1/20 + Y2 . (AlO) 



There is also a critical point (0,0), which corresponds to de Sitter spacetime with v = and = A/3 
if (|A2bp is used as a constraint equation. This critical point is not of interest to us now and further 
discussion of this case will be omitted. In addition, it can be shown that the critical point (0, 0) is not 
asymptotically stable. 
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Then, ( lASP becomes 
dYi 



-3 Fi + 3^2 = ^1, (Alia) 



ds 

dY2 yio + Yi 



A - {yio + Vl) 

X (a - (yio + Fi)' + 3 (yio + Fi) (^20 + F2) - 6 (1/20 + F2)') = G2 . (Allb) 

In order to prove that the critical points (?/io,l/2o) from (lASp are asymptotically stable 
solutions, it suffices to consider small Yi and Y2: \Yi\ \yio\ and IF2I ^ 1 2/20 1- We, then, 
find the following vector equation: 

^Y{s) = A-Y{s) + f{Y„Y2), (Al2a) 
ds 

with the vectors 

and the constant matrix 

1 / -9 +9 \ , ^ ^ 

The eigenvalues of A are both negative ((Xi = —2, cr2 = —4). The vector component /i is 
zero and /2 is quadratic in Yi or F2 to leading order: 

/i(Fi, r2) = 0, (Al4a) 
/2(Fi, 1^2) = O {Y,^ Yi, Fi Y2) , (Al4b) 
so that the following bounds hold: 

lim 4ffil^= lim 4E^ = o. (Al4c) 

With these results, the Poincare-Lyapunov theorem (Theorem 7.1 in Ref. [l^; see also The- 
orem 66.2 in Ref. [l5|) proves that the critical points (2/10,1/20) from (lASP are asymptotically 



stable (attractor) solutions. 



4. Stability analysis: Lyapunov function 



For completeness, we also give another proof which directly starts from ( lAlip . This proof 



relies on the construction of an appropriate Lyapunov function V , in order to be able to 
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apply the second Lyapunov stability theorem (l^ . Il5 |. The construction proceeds in three 
steps. 

The first step is to define the Lyapunov candidate function V[s, Yi, Y2] with properties 
V[s,0,0] = and V[s, Fi, ¥2] > for (Fi, Y2) ^ (0,0). Specifically, take the following 
quadratic function: 

V[s, Y^, Y2] = Y^ + Y^ + (Fi - ^2)' , (A15) 

which has the required properties and no explicit dependence on s. 
The second step is to calculate the orbital derivative, 

dV dV dV 

where the first derivative on the right-hand side vanishes for the choice ( lAlSp and where Gi 
and G2 are defined by the right-hand sides of (1A11I) . The explicit result for LgV is 

2 



L,V = - 



19 (Fi - Y2Y + 8 F2' + 24 (Fi^ Y2 - Fi^ + Fa^) 



4-(l + 2Fi)2 
-28 Fi^ + 56 Fi^ F2 - 60 Fi^ Fs^ + 48 Fi Fs^ 



with Taylor expansion 

2 
~3 



19(Fi-F2)' + 8F22 +0(f3). 



(A17) 



(A18) 



The third and last step is to demonstrate that (1A17I) implies the following inequality for a 
sufficiently small domain of Fi and F2 around Fi = F2 = 0: 



(n, 12)7^(0,0) 



< 0, 



(A19) 



where the strict inequality holds with the point Fi = F2 = excluded. Result (]A19|) implies 
that the function flAlSI) is a genuine Lyapunov function. 

The result (]A19|) for Lyapunov function (lAlSP now establishes the fact that the solution 



15|) 



F(s) = is asymptotically stable (Theorem 8.2 in Ref. [l^ and Theorem 25.2 in Ref. 

The precise mathematical definition of this asymptotic attractor behavior can be found 
in, e..g. Sec. 5.2 of Ref. [l4| (for a general discussion, see also Sec. 35 of Ref. 15|). In 



short, arbitrary starting values (fo(so), vq{sq)) in a sufficiently small domain give a solution 
(f(s), /i(s)) which asymptotically approaches the solution flA9l) as 'time' s runs towards 
infinity. 
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Appendix B: Attractor solution in a model with an Fa term 



ODEs and new variables 



In Ref. [7|, we considered a single- vector-field model with the combination [A + Co (Qs)"^ + 
kR'^A'^] in the action density for A > 0. For the case of Co = 1 and k, = —1/2, this 
corresponds to having a vacuum-energy-density term e = in the action (12 .lap , with 
Ba{x) = and Fa defined by (E^a]). 

The resulting dimensionless ODEs read: 

v + 3hv-Qh^v = 0, (Bla) 
= \ - (^ij + ShvY + 3hv{hv + 2v) +rM, (Bib) 
2h + 3h'^ = X- {v + 3hvY - [{Ah + 9 h^) - 2v^ - 2vv - Ahvv] -^tm, (Blc) 
rM + 4:hrM = 0, (Bid) 
for A > 0. With /i(r) = d(r)/a(r), the solution of ( IBldj) is known to be rjv/(T) oc l/a(r)^. 



The system of differential equations can then be rewritten as follows: 

V + 3hv -Qh^v = 0, (B2a) 



3{l + 2v^)h^-X + V^- TMstart (astart/«) = , (B2b) 

{l + 2v^)h~V^ + AhvV+- TAfstart (astart/o) ^ = , (B2c) 

where a = a(r) is the scale factor of the flat RW metric, agtart its value at r = Tgtart) and 
r Af start the corresponding starting value of tm- 
Next, introduce new variables: 

yi = V , 1/2 = hv, (B3a) 

s = ln(a) - In(astart) • (B3b) 



Writing ( ]B2p in terms of the new variables (]B3|) gives the following nonautonomous system 
of differential equations: 

y'l = Fi{s, yi, y2) , y'2 = F2{s, yi, 1/2) , (B4) 

where the prime stands for differentiation with respect to s and 

Fi{s, yi, 7/2) = -3 ivi - 2y2) , (B5a) 

ip. X , 3 2/1 - 12^/1^2 - 2 TAf start exp[ -4 s] mp-UN 

F2{s, yi, 2/2 = 2/1 + 2/2 — r~r~\ — • (^^^) 

A - 2/1 + rAf start exp[ -4 s \ 
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2. Critical points 



= lim ^2(5, yi, 1/2) 

y 10, 1/20 s^oo 



= 0. (B6) 

?;io, j/20 



A critical point (?/io,l/2o) of system (]B4p is defined as follows: 
lim Fi(s, yi, ^2 

s— ^00 

A straightforward calculation gives two such critical points, 

yt, = ±^/2X/5, yt, = ±^/X/10, (B7) 
corresponding to the asymptotic solutions 

Vty^^ = ±y^ T, /lasymp=(l/2)r-S (B8) 

in terms of the original variables.^ 

3. Stability analysis 

Make the following shift of variables: 

yi = yio + yi, y2 = y2o + Y2 . (B9) 

Then, (IB4p becomes 
dYi 



ds 

dY2 
ds 



-3 Fi + 6 F2 , (BlOa) 

yio + Yi + (1/20 + Y2) 

^ 3 {yio + Yif - 12 {yio + Yi){y2o + I2) - 2 rMstart exp[ -4 s ] 

A - {yio + + TAf start exp[ -4 S ] 

In order to prove that the critical points {yio, 2/20) from flB7p are asymptotically stable 
solutions, it suffices to consider small Yi and ¥2'. I^il ^ Iz/iol and IF2I ^ |l/2o|- We, then, 
find the following vector equation: 

Y{s) = A ■ Y{s) + B{s) ■ Y{s) + f{s, Y^, Y2) , (Blla) 

with the vectors 



^ It is obvious that (0, 0) is also a critical point, independent of the value of rMstart- See Ftn. [T]for further 
comments. 
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and the matrices 

l/-9 +18 \ ^, , 20a(s) , - - , ^ 

in terms of the auxihary variable 

a{s) = r A/start exp[ -4s] . (Blld) 

The eigenvalues of A are both negative (o"i = —4, (T2 = —5) and the matrix B{s) vanishes 
as s — )■ +00. The vector component /i is zero and /2 is quadratic in Yi or Y2 to leading 
order, 

Ms, Yu Y2) = 0, (B12a) 

. ( y,^ — i-. ^i^<i>±^ n r.) . (Bi2b) 

V(5a(s) + 3A) 5a(s) + 3A (5a(s) + 3A) / 

so that the following bounds hold: 

hm fl(^l^lM= hm M£l^lS = 0. (B12c) 

With these results, the Poincare-Lyapunov theorem (Theorem 7.1 in Ref. [l^; see also the 
discussion below Theorem 66.2 in Ref. 



15| ) proves that the critical points {yiQ, 1/20) from (IB 71 



are asymptotically stable (attractor) solutions. The above discussion provides the details 
for the result announced in the Note Added of Ref. [3| . 
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